The Chameleonic Contribution to the SZ Radial Profile of the Coma Cluster 
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We constrain the chameleonic Sunyaev-Zel'dovich (CSZ) effect in the Coma cluster from measure- 
ments of the Coma radial profile presented in the WMAP 7- year results. The CSZ effect arises from 
the interaction of a scalar (or pseudoscalar) particle with the cosmic microwave background in the 
magnetic field of galaxy clusters. We combine this radial profile data with SZ measurements towards 
the centre of the Coma cluster in different frequency bands, to find ATsz,rj(0) = —400 ± 40 fiK 
and ATcsz GHz (0) = -20 ± (68% CL) for the thermal SZ and CSZ effects in the cluster re- 

spectively. The central value leads to an estimate of the photon to scalar (or pseudoscalar) coupling 
strength of g c g ~ (5.2-23.8) x 10 _10 GeV _1 , while the 95% confidence bound is estimated to be 

3off 



« < (8.7-39.4) x HT 10 GeV" 



PACS numbers: 



I. INTRODUCTION 



There has been a great deal of interest over the past 25 
years or so in the mixing of axion-like particles (ALPs) 
with photons in the presence of magnetic fields; for exam- 
ple Refs. [l|-lH or see Ref. Q for a recent review. ALPs 
are a generalisation of the Peccei-Quinn axion originally 
introduced to solve the strong CP problem [(| . They re- 
fer to any light scalar or pseudoscalar with a linear cou- 
pling to the electromagnetic (EM) tensor, and have two 
free parameters: the axion mass and the EM coupling 
strength. For more details of the axion and the experi- 
mental bounds placed on its parameters see for example 
0, 0]- More exotic ALPs have been postulated, such 
as the chameleon scalar field developed by Khoury and 
Weltman Q. The chameleon scalar field has a density- 
dependent mass, making it an attractive candidate for 
dark energy. In high-density environments such as on 
Earth the chameleon is heavy and evades fifth-force lab- 
oratory searches despite a strong coupling to matter; in 
low-density environments the chameleon is light and can 
drive the cosmic acceleration. Bounds placed on the ax- 
ion in the laboratory are evaded by the chameleon and 
so the chameleon-photon coupling strength is less well 
constrained. The strongest bounds come from mixing 
in large-scale astrophysical magnetic fields, see for ex- 
ample Refs. [l], [t| UM- It was pointed out in pj that 
in sparse astrophysical environments, such as that found 
in galaxies or galaxy clusters, the mixing between the 
chameleon field and photons is indistinguishable from 
the mixing of any very light ALP. In our analysis we 
assume full generality as to the nature of the ALP but 
the resulting constraints will be most applicable to the 
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chameleon scalar field. A detailed discussion of the al- 
lowed chameleon parameters is given in Rcf. 11]. The 
best direct upper-bounds on the matter coupling 1/M 
come simply from the requirement that the chameleon 
field makes a negligible contribution to standard quan- 
tum amplitudes, roughly 1/M < Kr 4 GeV -1 [H El- 
Stronger constraints have been derived on the photon 
coupling g c g = 1/M c g. Laser-based laboratory searches 
for general ALPs, such as PVLAS [HI], and chameleons, 
such as GammeV pj|, give g e g < 10 _6 GeV _1 . Con- 
straints on the production of starlight polarization from 
the mixing of photons with chameleon-like particles in 
the galactic magnetic field provide an upper-bound of 
g eS < 9 x lO-^GeV" 1 Q. 

In [l(| we predicted the existence of a chameleonic 
Sunyaev-Zel'dovich (CSZ) effect arising from the mixing 
between a light scalar (or pseudoscalar) field and the cos- 
mic microwave background (CMB) radiation in the mag- 
netic field of galaxy clusters, with particular reference 
to the chameleon model. The scalar-photon interaction 
causes a fraction of the photons to be converted along 
the path, decreasing the overall photon intensity. The 
effect is similar to the thermal Sunyaev-Zel'dovich (SZ) 
effect arising from scattering of photons off electrons in 
the cluster atmosphere, but with a significantly different 
frequency dependence. In [l(| we compared predictions 
of the thermal SZ and CSZ effects in the Coma cluster of 
galaxies to the measured SZ decrement towards the cen- 
tre of the cluster in a number of frequency bands. The 
results constrained the scalar-photon coupling strength 
to be, geff < (0.72 - 22) x Kr 9 GeV~\ depending on 
the model that is assumed for the cluster magnetic field 
structure. 

In [l(| we proposed that the constraints would be 
stronger if the radial profile of the SZ effect could be anal- 
ysed. The thermal SZ effect scales with the integrated 
electron density in the cluster and decreases rapidly as 
the electron density decreases towards the edge of the 
cluster. The chameleonic SZ effect, by contrast, depends 
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on the ratio of the magnetic field to the electron den- 
sity, which introduces an extended wing to the predicted 
radial profile. This additional power in the SZ signal 
at larger radii is a distinctive signature. Until recently 
the only analysis has been of the stacked SZ profiles of 
galaxy clusters at known positions on the sky, suggest- 
ing a discrepancy between the SZ intensity decrement as 
measured by WMAP and the expected thermal SZ signal 
based on the X-ray emission data [HI, HH . This discrep- 
ancy takes the form of a greater than expected SZ signal 
at large radii. One possibility is that a simple beta-model 
is inadequate for describing the radial profile of the elec- 
tron density distribution in the cluster and that a more 
sophisticated model is required, but another more excit- 
ing possibility is the presence of a chameleon scalar field. 

The recent WMAP 7-year SZ measurements of the 
Coma cluster radial profile [l|| are the first detailed ob- 
servations of the SZ profile around a cluster to date. We 
calculate the CSZ contribution to the signal in the com- 
bined V + W band extracted from the WMAP measure- 
ments; we believe it is beyond the scope of this paper to 
perform a full analysis of the raw WMAP data contami- 
nated by noise and CMB anisotropies. 

This paper is organised as follows: in section |TT] we 
present the calculations for the mixing between ALPs 
and photons in the magnetic field of galaxy clusters. In 
section lTnl we discuss measurements of the magnetic field 
structure and SZ effect in the Coma cluster. The results 
of our analysis are in section IIVI and we conclude with a 
discussion of the results in section [Vj Appendix [A] con- 
tains the derivation for the probability of conversion be- 
tween photons and ALPs in galaxy clusters. Appendix 151 
outlines our simplistic approach to calculating the con- 
tribution of the CSZ effect to the V + W band in the 
WMAP data. 



II. MIXING BETWEEN PHOTONS AND 
AXION-LIKE PARTICLES IN CLUSTER 
MAGNETIC FIELDS 

The interaction Lagrangians describing coupling be- 
tween photons and scalar- or pseudoscalar-ALPs are, 

^scalar '^Qeffi&Fni/F 
-^pscudoscalar ^Qeffi&Fni/F : 

where F^ v = e^a^F 01 ^ /2 is the dual of the field strength 
tensor. The difference between the scalar and pseu- 
doscalar interaction determines which photon polariza- 
tion state is involved in mixing. The modification to the 
intensity of any radiation is therefore identical for cither 
particle but the modification to the polarization state is 
particle specific. 



The resulting field equations for this interaction are, 
n<p =m%4> +\g cS F^F»\ (1) 
(--^F^F^ , 

and, 

W V F^ =J» -g eS V v {ct>F^), (2) 

where the bracketed term is for a pseudoscalar as opposed 
to a scalar ALP. J M is the background electromagnetic 
4-current, such that V M J M = 0, and is the mass of 
the ALP. The Lagrangian for the chameleon scalar field 
must include the coupling between the scalar field and 
matter. However in sparse astrophysical environments 
the varying properties of the chameleon can be encapsu- 
lated in its mass to give the same field equations as 
Eqns. (JTJ and ©. See jlO] for a full derivation. 

We follow a similar procedure to that given in Refs. [l|- 
H, [l(| to determine the evolution of the ALP and photon 
fields, when propagating through a background magnetic 
field B(x). We expand the fields as perturbations about 
their background values, if = <f> — 4> and — A M — 
where F^ = d^A v — d v A^. Ignoring terms that are 
second order in the perturbations, we find 

-^ + V 2 ^m^+ 5cf[ B.(Vxa), (3) 
(+<7effB.(9 t a)) , 

and, 

- 9 t 2 a + V 2 a ~ w^a +g cS (V0 x B) , (4) 
(-. 9cff B(9 t 0)) , 

where again the bracketed term applies for pseudoscalars. 
The plasma frequency, w 2 j = 4iraEMn c /m c , arises from 
interactions between the photons and electrons of number 
density n c in the plasma, and acts as an effective photon- 
mass. 

The probability of conversion between photons and 
chameleons after passing throu gh t he magnetic field of 
a galaxy cluster was derived in [n|. The details are re- 
produced here in Appendix[Al and have been extended to 
apply generally to scalar- and pseudoscalar-ALPs. The 
derivation assumes weak-mixing between the photons 
and ALPs (valid for CMB photons propagating through 
galaxy clusters), and requires the ALP mass to be below 
the plasma frequency in the cluster. For the case of the 
Coma cluster this imposes < 10~ 12 eV. It was shown 
in [10( that for CMB photons travelling through the in- 
tracluster (IC) medium the mass of the chameleon scalar 
field satisfies this requirement. 

The magnetic fields of galaxy clusters have been mea- 
sured to have a regular component B rcg (x), with a typ- 
ical length scale of variation similar to the size of the 
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object, and a turbulent component <5B(x) which un- 
dergoes 0(1) variations and reversals on much smaller 
scales. Instead of considering a simple cell model for the 
magnetic field fluctuations, common in both the liter- 
ature concerning photon-ALP conversion in galaxy and 
cluster magnetic fields and in that concerning the mea- 
surement of those fields, we consider a more realistic 
model for B assuming a spectrum of fluctuations running 
from some very large scale (e.g. the scale of the galaxy 
or cluster) down to some very small scale (-C L co h)- 
We describe these fluctuations by a correlation function 
i?Bij(x) = (8Bi(y)SBj(x + y)}; here the angled brack- 
ets indicate the expectation of the quantity inside them. 
This description of the magnetic field is compatible with 
the requirement that V • B = [l7| • Assuming that the 
fluctuation spectrum is isotropic we have: 

RBij(x) = -5ijR^(x) + XiXjR^(x), 

where ii|(a;) is related to R%(x) by imposing V • B = 0. 
The photon-scalar conversion probability only depends 
on the components of the magnetic field that are per- 
pendicular to direction of the path, and so do not de- 
pend on We define Pe(fc) to be the power spec- 
trum associated with R^(x). We allow for a similar 
spectrum of fluctuations in the electron number den- 
sity, n c (z) — h c + 8n e (z), with auto-correllation function 
Rn{x) and associated power spectrum P^(k). 

The power spectra, Pn{k) and Ptss(k), have not been 
measured for the Coma cluster, but we draw on obser- 
vations of our own gala xy and other clusters to infer 
their probable form [1JJ, ll7Hl9l |. These measurements 
are generally consistent with a Kolmogorov-type turbu- 
lence model, but only probe the power spectrum at spa- 
tial scales larger than a few kiloparsecs. The dominant 
contribution to the probability of photon-ALP conver- 
sion in galaxy clusters is sensitive to Pe(fc) and P^(k) 
on scales < 10 -3 -0.1pc. There are no observations 
of the form or magnitude of Pb(&) and P^(k) on such 
small scales. According to Kolmogorov's 1941 theory 
of turbulence, the power spectrum of three dimensional 
turbulence on the smallest scales is universal and has 
P(k) oc fc~ n / 3 . We therefore estimate Pb(A:) and P^(k) 
in the Coma cluster on small scales by assuming we have 
three dimensional Kolmogorov turbulence. 

Wc find that the dominant contribution to the average 
probability of photon-ALP conversion in a galaxy cluster 
is (Eq. ([Al3]) . Appendix El, 

(5) 

where 1-^ = 1 + (5ri^/n^), v is the frequency of the ra- 
diation, L is the path length through the magnetic re- 
gion, and Lb is the coherence length of the magnetic field 



which is defined later in Eq. The range depends on 
the normalisation for the power spectra Pe(fc) and P^(k). 
The resulting modification to the CMB intensity is, 

AI _ - 
Jo 

which holds equally for both scalar- and pseudoscalar- 
ALPs. 



III. THE COMA CLUSTER 

The nearby Coma cluster of galaxies is unique in hav- 
ing detailed measurements of both its magnetic field 
structure and the SZ effect. In this section we discuss 
these measurements. 



A. Magnetic field 

The magnetic field at the centre of the Coma clus- 
ter was determined by Feretti et al. [20] from observa- 
tions of Faraday rotation measures (RMs) of an embed- 
ded source. The RM corresponds to a rotation of the po- 
larization angle of an intrinsically polarised radio source. 
The RM along a line of sight in the z direction of a source 
at z — z s is given by: 

f Zs 

RM(z s z) = ao / n a (xz)B\\ (xz)dx, 
Jo 

where ao = a^j^/n 1 ^ 2 !!! 2 , B\\ is the magnetic field along 
the line of sight, and the observer is at x = 0. For a 
uniform magnetic field over a length L along the line 
of sight and in a region of constant density, we expect 
(RM) = aoB»n e L. However, the dispersion of Faraday 
RMs from extended radio sources suggests that this is 
not a realistic model for IC magnetic fields. Both the 
dispersion of the measured RMs and simulations suggest 
that IC magnetic fields typically have a sizable compo- 
nent that is tangled on scales smaller than the cluster. 
The simplest model for a tangled magnetic field is the 
cell model where the cluster is divided up into uniform 
cells with sides of length L co h- The total magnetic field 
strength is assumed to be the same in every cell, but the 
orientation of the magnetic field in each cell is random. 
In the cell model the RM is given by a Gaussian distri- 
bution with zero mean, and variance 

<RMV%^V(z))n c 2 (z)dz, 

where (B 2 ) is the average value of B 2 along the line of 
sight. A more realistic model for the magnetic field is to 
assume a power spectrum of fluctuations. In this model 
one finds fl9j ]. 

(RM 2 )^^^(B 2 (z))n 2 (,)dz, 
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where following Ref. [19j the correlation length scales, 
Lb and Ln, for the magnetic and electron density fluc- 
tuations respectively are defined, 



J™kdkP B/N (k) 

2/ oo fc2 dA ,p B/N(fc) - 



(6) 



To match the same observations with the same total field 
strength, Lb in the power spectrum model should be 
equal to L co h in the cell model. 

Assuming a cell model for the Coma cluster magnetic 
field in combination with a uniform magnetic field, Fer- 
etti et. al. 20] detected a uniform component of strength 
0.2 ± 0.1 /LiG coherent on scales of 200 kpc and a tangled 
component of strength 8.5 ± 1.5 /xG coherent on much 
shorter scales of 1 kpc extending across the core region of 
the cluster. 



B. Radial Profile 

The electron density distribution plays a crucial role 
in determining both the radial profile of the thermal SZ 
effect and that of the CSZ effect. The electron density of 
cluster atmospheres can be determined from X-ray sur- 
face brightness observations, for example using ROSAT, 
which are generally fitted to a simple beta profile (see 
Ref. [HI): 



n e (r) = n a 



r 2\ -■> ; - 



(7) 



where hq is the electron density at the cluster centre, r c 
is the core radius (a free parameter fitted to the X-ray 
observations) and j3 ~ 0.6-1 in general. For the Coma 
cluster the best fitting parameters are (22|, 

P = 0.75 ±0.03, 
6 c =r c /Dp L = 10.5 ± 0.6 arcmin, 

where Da = 67 h^ 1 Mpc is the angular diameter distance 
to the Coma cluster. The average electron density at the 
centre of the cluster has been measured at 4 x 10 -3 cm~ 3 . 

In a similar fashion the magnetic field is expected 
to decline with radius and the coherence length to in- 
crease. Magneto-hydrodynamic (MHD) simulations of 

1/2 

galaxy cluster formation suggest that (B 2 ) oc (n e ) ri 
with 77 = 1 while observations of Abell 119 suggest 
7] = 0.9 [23j]. Two simple theoretical arguments for the 
value of rj, are (1) that the magnetic field is frozen-in 
and hence r\ = 2/3, or (2) the total magnetic energy, 
oc (B 2 ), scales as the thermal energy, oc n e , and hence 
i] = 1/2 [2l|. Faraday RM observations of clusters [24| 
have found that the magnetic field extends as far as the 
ROSAT-detectable X-ray emission (coming purely from 
the electron density distribution). Simulations have also 
shown that the coherence length increases with distance 
from the cluster core, however the radial scaling of L C oh 



Experiment 


^/GHz 


Sp/GBz 


ATsz/fiK 


OVRO 


32.0 


6.5 


-520 ± 83 


WMAP 


60.8 


13.0 


-240 ± 180 


WMAP 


93.5 


19 


-340 ± 180 


MITO 


143 


30 


-184 ±39 


MITO 


214 


30 


-32 ± 79 


MITO 


272 


32 


172 ± 36 



TABLE I: SZ measurements of the Coma cluster core in dif- 
ferent frequency bands. 



is uncertain. We make the simple ansatz that, at least 



inside the virial radius, L co h oc n, 
that 



-7/3/3 



for some 7, so 



L, 



coh 



L 



cohO 



r 2 \ 7 -> 



where L C oho is the coherence length in the centre of the 
cluster. Based solely on dimensional grounds, estimated 

— 1/3 

values for 7 are: 7 = p so that L co h oc n c , or 7 = 1 so 
that L C oh oc r for large r. The simulations of Ref. [25| . 
suggest that for large r, L co h roughly scales as r inside 
the virial radius. 



C. Sunyaev-Zel'dovich measurements 

The thermal SZ effect arises from inverse- Compton 
scattering of CMB photons with free electrons in a hot 
plasma. The scattering process redistributes the energy 
of the photons causing a distortion in the frequency spec- 
trum of the intensity. This is observed as a change in the 
temperature of the CMB radiation in the direction of 
galaxy clusters [III [2(| , 



77l fi 



T, 



11 



(8) 



where the integral is along the line of sight. The fre- 
quency dependence of the effect is contained in gsz {v) = 
xcoth (I) — 4 where x = hv/ksTo. At low frequencies 
the SZ effect is seen as a decrement in the expected tem- 
perature from the CMB, while at high frequencies it is 
seen as a boost in the temperature. The electron pres- 
sure, P c = n c k,B,T c , varies with the distribution in the 
cluster, where T e is the temperature of the electrons in 
the plasma. The SZ effect is also commonly expressed in 
terms of the optical depth, defined by To = J aTn c (l)dl. 

Measurements of the CMB temperature decrement to- 
wards the centre of the Coma cluster were assimilated 
and analysed in (27j . These measurements are summa- 
rized in Table |TJ Assuming only a thermal SZ contribu- 
tion the optical depth was found to be 4.7 x 10~ 3 . This 
assumed an electron temperature of k^Tc = 8.2 keV. 

An alternative parameterisation of the thermal SZ ef- 
fect is to consider the central temperature decrement in 
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the Rayleigh- Jeans (RJ) limit (y — y 0). Separating out 
the frequency and radial dependence, and assuming a 
beta model for the electron distribution (Eq. (J7J)), we 
have 



AT SZ {9) = AT sz , RJ (0) 



°1 
9 2 



1 (l"3/3)/2 



where the radius r from the centre of the cluster to some 
point in the cluster atmosphere has been approximated 
by r 2 ks (DaO) 2 + I 2 . Here 9 is the angle subtended 
from the cluster centre to the point and I is the radial 
distance of the point from the cluster centre along the 
line of sight. Within this parameterisation the inferred 
optical depth for the Coma cluster, tq w 4.7 x 10~ 3 , 
corresponds to a central temperature decrement in the 
RJ limit of -420 fiK. 

Detailed measurements of the CMB temperature in 
43 radial bins centered on the Coma cluster have re- 
cently been released by the WMAP collaboration [l6j . 
Assuming a beta model for the electron density distribu- 
tion, they find a central temperature decrement in the 
Rayleigh- Jeans limit of —381 ± 126 fiK from the data 
extracted in the V band (60 GHz), and -523 ± 127 
from the W band (90 GHz). This discrepancy is most 
likely caused by contamination from the primary CMB 
anisotropies. The data is processed to eliminate this 
CMB signal by combining the V and W bands, to give 
AT S z,rj(0) = '-377± 105/iK (68% CL). This suggests a 
CMB anisotropy signal of the order — 100 /zK coinciding 
with the centre of the Coma cluster. The extracted data 
in the combined V + W band is reproduced in Fig. [1] 



V+W combined 




100 150 200 250 300 

8 / arcmin 



FIG. 1: Measurements of the cleaned SZ signal for the Coma 
cluster in the combined V + W band, plotted in units of the 
Rayleigh-Jeans (RJ) temperature, courtesy of Eiichiro Ko- 
matsu [l(|. The solid black line is the best fit thermal SZ 
model assuming a beta profile for the electron distribution 
with no CSZ contribution. The dashed blue line is the best 
fit model assuming both thermal SZ and CSZ contributions. 



IV. THE CHAMELEONIC 
SUNYAEV-ZEL'DOVICH EFFECT IN THE 
COMA CLUSTER 

We saw in section [XT] that the average conversion rate 
of photons into ALPs after passing through the magnetic 
field of a galaxy cluster was given by Eq. (JQ). For the 
Coma cluster the regular component of the magnetic field 
has been measured to be much smaller than the tangled 
component (true of most cluster magnetic fields) , and so 
the conversion probability is approximately, 



V. 



1^4 



i. ( „,™. 45)& (^ 



-5/3 



(SB 2 ) I 2 (21* 



1 



A reasonable value for In = 1 + (5n 2 /n 2 ) based on elec- 
tron fluctuations in our own galaxy is In ~ 1—2. The 
photon- ALP conversion causes a decrement in the CMB 
temperature towards galaxy clusters, with 



AT, 



CSZ 



- 1 



To 



-R 



where x = hv/k-oT^. The above equation assumes the 
average magnetic field, electron density and coherence 
length are constant along the path length L. However 
if the fields vary slowly over the path length, then the 
effect will be proportional to the integral along the line 
of sight: 

P W oc J dl(SB 2 )n- 5 / 3 L B 2/3 . 

Assuming a beta model for the electron density in the 
cluster and the radial scaling for the magnetic field and 
coherence length discussed in section UlI Bl we have 



AT CS z(0) 



A rp20i GHz 
^ J CSZ 



1 



(0) 



x e~ 



5/3 



IyiO) 



(9) 



where v = 204 GHz and we have defined p = 1 — 27/3 - 
677/3 + 5/3, and 



(1 



+ y 2 ) 2 w(y)dy, 



with y 2 = (l/r c ) 2 ^1 + j^j . We assume the turbulent 

cluster field dies off very quickly beyond some charac- 
teristic radius r max , due to Lb growing very quickly or 
((5B 2 (z)) decreasing much faster than n^. A reasonable 
assumption for r max is to approximate it by the virial 
radius r v - lr [10( . Thus we have introduced a window func- 
tion, w = exp(— r 2 /r^ ir ), which corresponds to 

+ (gg + 9 2 )y 2 



w(y) = exp - 



6 



where 6 V [ Y = t v -„IDa- The magnitude of the CSZ effect 
at the centre of the cluster at 204 GHz is given by, 



ATA 



z (o) 



- (0.13-2.69) x W- 7 gl 

-5/3 



SB 



2\l/2 



8.5 fj,G 
2/3 



(ixTIF^cnF 7 ) (l50kpc) (lkpc) ' (^^) 

where g w = W(10~ 10 GeV- 1 ). 

The WMAP measurements of the Coma cluster radial 
profile discussed in section IIII CI constrain the degree of 
photon-ALP mixing in the Coma cluster. We calculate 
the contribution of the CSZ effect to the combined V + W 
band data. This assumes a simplistic approach to the 
processing technique of Komatsu et al. [16| in combining 
the V and W bands to eliminate contamination from the 
CMB anisotropics. This is detailed in Appendix [B] We 
predict the total signal to be, 



(0) 



AT. 



SZ.RJ 



(0) 



. a t-i204 GHz 



(0) 



2 (.9 csz 



3csz) 



5csz(fo)(5sz 



9&) 



where gcsz{v) 



.5/3 



(e- x - l)/x. 



Comparing our predictions, T^ cd , to the extracted 

WMAP data in the combined V + W band, 
the 43 radial bins of the Coma cluster, we maximise the 
likelihood L defined by, 



T ex P t> over 



-2\ogL = J2 



T. 



cxpt 



prod 



cxpt 



where er^ pt are the errors in the WMAP measurements. 
We take AT sz ,rj(0) and AT^| GHz (0) to be the free pa- 
rameters in the theory and assume no prior knowledge 
of their distributions. We marginalise over the remain- 
ing model parameters 6 C , 7, rj, v ir), and assume the 
following priors on their distributions: 

13 ~ N (0.75, 0.03), 

9 C ~ A/" (10.5, 0.6) arcmin, 

7 ~ U (0.6,1), 

7 1 ~ W (0.5,1), 

6» vir - A/" (102, 20) arcmin, 

where A/" (/1, cr) is a normal distribution with mean \x and 
standard deviation a; and W (a, &) is a uniform distribu- 
tion over the range a to b. These priors come from the 
experimental measurements and theoretically motivated 
values discussed in section ITlI B I The virial radius of the 
Coma cluster is 1.99 ± 0.22ft, _1 Mpc [28[ which, coupled 
with Da ~ 67ft.~ 1 Mpc, gives the above prior on 9 v - lr . To 
find the confidence limits on our results we assume that 
X 2 = —2 log (li/Lj follows a xi distribution, where L is 
the maximised likelihood over the parameter space. 



The results of this maximum likelihood procedure are 
presented in Fig. The blue dashed line shows the 68%, 
95% and 99% confidence limits on the parameter space 
of AT sz ,rj(0) and AT2°| GHz (0) . Considering only these 
WMAP radial measurements we then find that, 



AW m (0) = -60±25 M K 



n204GHz/ 
CSZ \ 

AT sz ,rj(0) = -480 ± 100 nK (68%CL). 
In addition to this data set, the SZ measurements to- 
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FIG. 2: Confidence limits on the parameter space of 
AT sz ,Rj(0) and AT ( §t GHz (0) for the Coma cluster. The 
dashed blue lines are the 68%, 95% and 99.9% confidence 
limits obtained from the WMAP 7-year measurements of the 
SZ radial profile in the V + W band. The red dotted lines 
are the confidence limits obtained from the SZ measurements 
towards the centre of the cluster in different frequency bands 
first presented in 10] . The black solid lines are the combined 
confidence limits from these two data sets. 

wards the centre of the Coma cluster given in Table U 
were analysed by the authors in [10( . The resulting con- 
straints on the parameter space were, 



ATa° s 4 z GHz (0) 
AT sz ,rj(0) 



> -20^iK 

= -420±60/;K (68%CL), 



and are reproduced in Fig. [5] as the red dotted line set. 
Combining these two data sets we find, 

AT 204GHz( ) = _20±15AiK 



CSZ 

AT sz ,rj(0) 



= -400±40^<K (68%CL), 



where the 68%, 95% and 99.9% confidence limits on the 
parameter space are also plotted in Fig. [2] We note 
that, assuming no CSZ contribution, Komatsu et al. [HJ 
inferred AT sz ,rj(0) = -377 ± 105 fiK. 

The CSZ contribution to the WMAP V + W band, 
implied by AT^f z Giiz (0) = -20 fiK, is plotted in Fig. ffl 
In Fig. [3] we have plotted the expected contributions to 
the signal at 60 GHz, 90 GHz and 214 GHz. The thermal 
SZ effect is approximately null at 214 GHz and so future 
measurements at this frequency will be most affected by 
the CSZ contribution. 
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FIG. 3: The predicted contributions from the CSZ (solid red line) and thermal SZ (dashed black line) effects in the Coma 
cluster at three different frequencies. We have taken the central temperature decrements to be those of the maximum likelihood 
parameters, AT sz ,rj(0) = -400rtK and ATcg| GHz = -20^K, and have assumed 6 C « 10.5arcmin, 7 = /3 w 0.75, rj m 0.9 and 
Ovii ~ 102 arcmin. 



The maximum likelihood estimate, ATgg| GHz (0) = 
—20 /iK, corresponds to a photon-scalar conversion prob- 
ability in the Coma cluster of "P 74 ^(204 GHz) = 2.7 x 
1CP 5 . Substituting into Eq. (fit)]) , we find an implied cou- 
pling strength between photons and axion-like particles 
of, 

g cS w (5.2-23.8) x HT 10 GeV _1 . 

Similarly the 95% confidence bound on the CSZ con- 
tribution from the combined data set, AT^\ Giiz (0) > 
—55 fiK, corresponds to ^ 7 ^(204 GHz) < 7.5 x 1(T 5 , 
giving a bound on the photon-scalar coupling strength 
of, 

9cs < (8.7-39.4) x KT^GeV^ 1 (95% CL). 

V. DISCUSSION 

In this work we have constrained the degree of conver- 
sion between photons and ALPs in the magnetic field of 
the Coma cluster. The existence of the ALP induces a 
chameleonic SZ effect which was first predicted in Ref. 
[To| . Measurements of the SZ effect in the direction 
of the Coma cluster have been made by a number of 
CMB probes. Considering only the measured radial pro- 
file of the SZ effect in the combined V + W band, we 
find AT£° s 4 z Giiz (0) = -60 ± 25^K, for the change to the 
CMB temperature due to the CSZ effect towards the cen- 
tre of the cluster. Considering instead SZ measurements 
towards the centre of the cluster in different frequency 
bands, the constraint is AT<?g| GHz (0) > -45 at 95% 
confidence. Combining these results we find, 

AT 204GHz( ) w _20±15^K (68% CL), 
and AT^| GHz (0) > -55 fiK (95% CL). 

The above result involves a number of assumptions. 
Most importantly, it assumes that the electron density 
distribution in the cluster can be modelled by a simple 
beta profile (Eq. [7]), common in both the measurements 
of cluster surface brightness and SZ profiles. However 



measurements of the stacked SZ profiles around clusters 
have suggested that this is a poor fit to the data and that 
a more sophisticated model is required [HI, [l6| . 

As was pointed out in [loj , estimates of the correspond- 
ing photon to scalar (or pseudoscalar) coupling strength 
is highly dependent on the model that is assumed for 
the cluster magnetic field. In this analysis we have al- 
lowed for a spectrum of fluctuations in the magnetic field 
and electron density and assumed that fluctuations on 
small scales can be approximated by three-dimensional 
Kolmogorov turbulence. This leads to an estimate of the 
photon to scalar coupling strength for the maximum like- 
lihood scenario of, 

g cS ks (5.2-23.8) x 10- 10 GeV~\ 

while for the 95% confidence bound on the coupling 
strength is, 

g c s < (8.7-39.4) x Hr 10 GeV -1 . 

The range depends on the magnitude of electron density 
fluctuations in the cluster. 

In our analysis we have also imposed that the ALP 
mass be less than 10 _12 eV. Combining this with the 
above coupling strength, we find that a standard ALP 
with those properties is ruled out, since it would violate 
constraints on ALP production in the Sun, He burning 
stars and SN1987A. However, these constraints feature 
ALP production in high density regions, and therefore 
do not apply to a chameleon scalar field. This includes 
both the chameleon discussed earlier and any other non- 
standard ALP for which one or both of g c s and 1 / de- 
crease strongly enough as the ambient density increases. 
The above estimate and bound on the chameleon-photon 
coupling strength is consistent with constraints that have 
been derived elsewhere [l|, H, [T(| • 

In our radial analysis of the Coma cluster we calculated 
the CSZ contribution to the V + W band extracted from 
the WMAP data. To confirm our findings, future work 
should include an analysis of the raw CMB data in the 
separate V and W bands. As we can see from Fig. [3l the 
CSZ effect is predicted to dominate over the thermal SZ 
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effect at 214 GHz and so measurements at these higher 
frequencies will be most sensitive to the existence of a 
chameleon. 

We conclude by noting that the anomalous dip in the 
WMAP measurements at large radii that can be seen in 
Fig. [T] cannot be attributed to photon-scalar mixing: ex- 
tending the window function that imposed rapid decay of 
the CSZ contribution beyond the virial radius, merely de- 
creases the total SZ contribution uniformly and provides 
a poorer fit to the data. 
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ichiro Komatsu for helpful comments on a draft version 
of this manuscript, and are grateful to Eiichiro Komatsu 
for providing the SZ radial measurements of the Coma 
cluster in the combined V+W band. 



Appendix A: ALP-Photon Conversion in Cluster 
Magnetic Fields 

In this section we present the calculations for the prob- 
ability of conversion between photons and ALPs in the 
magnetic field of galaxy clusters given in Ec< (JS|). The 
bulk of this derivation was first presented in [1JJ] . We as- 
sume a Kolmogorov-type turbulence model for the mag- 
netic field structure in the cluster and allow for similar 
fluctuations in the electron density. 

Starting from the mixing equations for scalars and 
pseudoscalars, Eqns. and (g]) given in section HH we 
consider a photon field propagating in the z direction, so 
that in an orthonormal Cartesian basis (x, y, z) we have 
a = ("f x , j y , 0) T . The equations of motion for the scalar- 
ALP and photon polarization states can then be written 
in matrix form: 





( 













d 2 t +dl- 





-I 












\ M e ff z 


A/cff ° z 











addition to this, we define the state vector 

% (z) 

3 2lA ^(p{z) 



where 



A(z) = f 
Jo 



dx, 



and m^ ff = m^z) — Wpj(z), which simplifies the above 
mixing field equations for a and (p to 



u 



B(z) 



2M 



u. 



(Al) 



where 



B{z) 



( 

B x e 

\B y e 2lA ~B x e 2lA 



-B y e- 2lA 

-2iA 



in the case of scalar-photon mixing, and 

B r e- 2lA 



B(z) 



-B x e 2lA -B y e 2iA 



B y e- 2tA 




in the case of pseudoscalar-photon mixing. 

To solve this system of equations we expand, u = Uq 
Ui + U2 + . . ., such that 



B[z) 



2M 



cir 



Ui ; u' a = 0. 



and neglect the higher order terms from mixing. We 
define 



Mi 



B(x) 



dx. 



while for a pseudoscalar-ALP the mixing matrix is 



o 2M ff 

We say that mixing is weak when 
Q ix(M\(z)Mi(zj)<l. 

Thus Eq. (jAll) can be solved in the limit of weak-mixing: 

u{z) ~ [I + M 1 (z) + M 2 (z)]u(0), 
where A4 2 = Jq M'i(x)A4i(x)dx, which can be written 



-d? + d 2 - 







o 



By 

M, 



ft 



■d. 



Afoff 

By 

' Meff ' 
™2 




= 0, 



u(z) 



I + Mi(z) + ±(M c {z) + Ml(z)) 



u(0), 



where we have written g c g — l/M c ff. Following a 
similar procedure to that in Ref. [3|, we assume the 
fields vary slowly over time and that the refractive index 
is close to unity which requires m^/2w 2 , uj 2 1 /2uj 2 and 
\B\/2ujM cS all < 1. Defining 7l = ^(z) e M*-i)+«>M 
and (p — tpe tu '^ z ~ t ' >+1 ^^ where /3 jZ = — lj 2 1 (z)/2lj, we ap- 
proximate — d 2 w id 2 and u> 2 + d 2 « 2oj(lj + id z ). In 



where explicitly in the case of scalar-photon mixing 

-Al 

Mi = 



M 



c 




a 



yy 
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and 



A, 



Cij — 



, BAx)e 2lA{ -^ 

ax . 

2M cff 



(A2) 



Ax (A*'(x)A j (x)-AUx)A' j (x)) 



The same holds for pseudoscalar-photon mixing with the 
replacements A x — > Ay, A y — > —A x , C xx O- C yy and 

C x y y Cy x . 

The polarization state of radiation is described by 
its Stokes parameters: intensity I 7 (z) = ^(z)! 2 + 
|7 y (z)| 2 , linear polarization Q(z) = \l x {z)\ 2 - \l y {z)\ 2 
and U{z) — 2Re (n x { z )ly{ z )) 1 an d circular polarization 
V(z) = 21m (7*(z)7 y (z)). Assuming there is no initial 
chameleon flux, 1^ = \(p\ 2 = 0, then to leading order the 
final photon intensity is given by 

L ( (z) = J 7 (0)(1-7W*)) 

+Q(0)Q q ( z ) + U(0)Q u {z) + V(0)Q v (z), 

where we have defined, 

*W*) = \{\A x {z)\ 2 + \A y (z)\ 2 ), 

<2 q (z) = ±(\A x (z)\ 2 -\A y (z)\ 2 ), 

Q u (z) = Re(A*A y ), 
Q v (z) = lm(A* x A y ). 

For the CSZ effect, we are concerned with CMB radi- 
ation, the intrinsic polarization of which is small, i.e. 
v/Q 2 (0) + U 2 (0) + V 2 (0)/I 7 {0) < 1. It follows that 
to leading order the modification to the intensity of 
the CMB radiation, from the effect of both scalar- and 
pseudoscalar-photon mixing, is 



AJ 7 (z) 



-7V^(z)/ 7 (0). 



We now proceed to evaluate the Ai assuming a 
Kolmogorov-type turbulence model for the magnetic field 
and electron density fluctuations. 

The magnetic field is split into a regular compo- 
nent, coherent on the scale of the cluster, and a fluc- 
tuating component which undergoes field reversals on 
a much smaller scale, B(z) = B rcg + (5B(z). We as- 
sume that, for each i and for fixed x, the turbulent 
component of the magnetic field SBi (x) is approximately 
a Gaussian random variable. We require (at least ap- 
proximate) position independence for the fluctuations 
so that i?Bij(x;y) sa R Bij (x), where i? B jj(x;y) = 
(SBi(y)SBj(x+y)). i?By(x) is the auto-correlation func- 
tion for the magnetic field fluctuations. 



We require that oj 2 { 3> m^, so that m 2 s 



V 



— n c . Spatial variations in m 2 s are then entirely due to 
spatial variations in n c . The electron number density 
is divided into a constant part and a fluctuating part, 
n c = n c + 5n c , and we assume that 1 + Sn e /n e is a log- 
normally distributed random variable with mean 1 and 



variance (S 2 ), where <5 n = 5n e /n e . We define the electron 
density auto-correlation function by, 

i? N (x) = (5n c (y)6n c (x + y)) . 

We assume isotropy, which implies J?n(x) = Rn(x). Sim- 
ilarly, isotropy allows us to write the magnetic field auto- 
correlation function in terms of two scalar correlation 
functions: 

fl B y(x) = ^R^(x)Sij + x i XjR B (x), 

where x is the unit vector in the direction of x. 
For example, taking x in the z-direction gives x = 

(0,0,1). We note that Tr (R Bij (x)) = R^(x) + R B (x) = 
(<5B(y) • (5B(x + y)}. The magnetic field must obey 
V • B = which implies 9ii?Br/(x) = and hence, 



3 dec x 



dx 



0, 



and we must have R B (0) = 0. This fixes R B (x) in terms 
oiR^(x). 

The power spectra for the magnetic and electron den- 
sity fluctuations, Pb(&) and P^(k), are defined: 

4W = ^ / d 3 fce 2 " k x P B (fc) 

= J fc 2 dfcP B (fc)sinc(27rfcx), 

R n (x) = J- / d 3 fce 27rlkx P N (fc) 

47T J 

k 2 dk P^(k)sinc(2irkx) , 



where sinc(a:) = sinx/x. 
We also define, 

I f z 

$n(z) = ~ S n (x)dx. 

z Jo 

Note S n (L) = 0,where L is the path length through the 
magnetic region. With Ai defined by Eq. (|A2[) , we define 



Z = (1 + S n (z))z and A = m 2 ff i/4w oc -n^L/Au). It 
follows that A(z) = AZ/L. We assume that \S n (z)\ < I 
along the majority of the path, and so Bi(z) pa Bi(Z). 
We then have, 



ML) 



dZ- 



B,{Zi) 



i ■ .a : 

-e i . 



, 2Af cff (l + <5 n (Zz)) 
We define the correlation, 

Aj = (Ai(L)^(L)) 
= ilfc Jo dx Io d y R \f d x ~ y)z)e 2 ^-y)/ L , 
where the correlation function is defined, 

Biiy^Bjdy + x)z) \ 



R^(xz) = 



(l + S n (yz))(l + 6 n ((y + x)z)) i 
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We split the magnetic field into a regular and turbulent 
component: B = B reg + <5B. Thus, 



This integral may be evaluated explicitly in terms of the 
function, 



R%\xz) = BT s Bp 

■<(l + « n ( W z))- 1 (l + 5 n (( 1/ + a; )S))- 1 ) 
SBi(yz) ■ SBj((y + x)z) 
(l + S n (yz))(l + d n ((y + x)z)) 

For simplicity, we assume that fluctuations in the mag- 
netic field and electron density are uncorrelated. Thus, 

■<(l + ( 5 n (y))- 1 (l + ^„(y + x))- 1 ). 

Defining, 

R S ( X ) = ((l + ,5 n (y))- 1 (l + 5 n ( y + x))- 1 ), 
and assuming isotropy we have, 

R\f(xz) = (V s Sf 5 + ^RHx) + ZiZjRl(x)j R 5 (x), 

where z = (0,0, 1). The photon-scalar conversion prob- 
ability is given by T 1 ^ (j} = \ (\A X \ 2 + \A y \ 2 ) = ±(z x 

A)t(z x A) where A = (A X (L), A y (L), A Z {L)) T . Thus 
the average conversion probability depends only on the 
components of A, and hence i?*°'(a;z), that are perpen- 
dicular to z. We therefore define, 

R\f (xz) = B^Bf s R 5 (x) + ^SijRkWRsix), 



and take, 



A 



" 4M c 2 ffJ0 , 



Since R^ B (x) must vanish when x — 0, we have 
Eff (0) = i?*f (0) = (bT s B^ + S f(\5B(yz)\^ R s (0). 
We define a combined power spectrum P*? t (fc) by, 



dx I dyRl?((x-y)z)e 



2iA(x-y)/L 



R%f(xz) = I k 2 dkPl°\k)smc(2nkx). 



We then have, 



r2 poo 

A — ' 1:1 r "" ( 



am 2 I ~ y 

4iW cff JO 



k 2 Pl°\k)K(TTkL- 1 A)dk, 



where, 



, r , i , . sin(27rfcLs) . T . 

KtnkL; A) = ds 2(1 - s) — ^ -cos(2As). 

' znkLs 



Si(x) 



sin y 



dy. 



/o y 

We define 5Si(aj, y) = Si(a;+y) — Si(a; — y). We then have, 

2xK(x;A) = (5Si(2A,2a;) 

A sin 2A sin 2x + x cos 2 A cos 2x — x 



x 2 - A 2 

We are concerned with CMB radiation, for which to ~ 
10 _5 -10~ 3 eV, propagating over distances of the order of 
100 kpc through a galaxy cluster. It is easy to see that 
for such a scenario |A| ^> 1. We evaluate K(x; A) in the 
asymptotic limit where |A| 3> 1. When x <C |A| we have, 



AxK{x; A) 



2x 
A 2 



and when x S> I A I , we have 



4xK(x; A) ~ 2vr 



sin 2a; cos 2 A 
A 2 ' 

2 

a; 



We define fc cr it = \A\/irL. We assume that the fluc- 
tuations in B and n c are such that P^° t (k) drops off 
faster than fc~ 3 for all k > fc* where fe* <C fc C rit- This 
is equivalent to assuming that the dominant contribu- 
tion to i?*°'(0) comes from spatial scales than are much 
larger than k~ rit . We expect that this dominant contribu- 
tion will come from scales of the order of the coherence 
lengths of the magnetic field and electron density fluc- 
tuations (Lb and Ln defined in Eq. ([B])) and so we are 
assuming that k~ lit <C Lb, Ln- With this assumption, in 
the limit of large |A|, we have to leading order, 

t2 poo 

±M 2 s Aj * ^ I k 2 P'°\k)dk 



2A 2 Jo 
L cos 2A 
4ttA 2 



/ fcsin(27r/fcL)i^ ot (fc)dfc 
Jo 



We define, 



t r°° 
+- kPt°\k)dk. 

qP^(q)dq, 



and we can then rewrite Aij as, 
A.. ~. 2uj2 r>'"< 



^cff™eff 



Ri°\0) 



2lo 2 



M 2 rfi 4 y 



Rlf(Lz) cos 2 A 



8M 



2~ Wjj (fccrit; 
off 



We have that, 



^f(o) = i?*f(o) « ( sr s ^ rcs + ^ (sb 2 ) e 
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We assume that fluctuations in n e follow an approxi- 
mately log-normal distribution, and hence 



r 3 



We assume, as is the case in the situations we consider, 
that L is much larger than the correlation length of the 
magnetic and electron density fluctuations. This implies 
that the only contribution to P*°*(Lz) which may result 
in a leading order contribution to Aij is, 

Finally, we must evaluate Wij(k) in terms of the mag- 
netic and electron density power spectra. We separate 
the fluctuations in n c into short and long wavelength 
fluctuations, S s and Si respectively, and assume they are 
approximately independent. Thus n c = n c (l + Si)(l+S s ). 
This parameterisation is consistent with the assumption 
that n e /n e has an approximately log- normal distribution. 
We assume that the short wavelength fluctuations are lin- 
ear up to some cut-off scale k^ . Above this spatial scale 
we have the long wavelength fluctuations, which are not 
necessarily linear. We then have that over small spatial 
scales, < jfej^, 



^(x) « <^^((l + <5 s (y))- 1 (l + 5 s (y + x))- 1 ), 

J)(f)~V^N(x)], 

= 4[l+n e - 2 # N (x)] . 
It follows that for k 3> fcnn we have approximately, 

P s (k) « I*n- 2 P N (k). 
We assume that fcn n <C fe C rit- Finally we define, 

P BS (k) = P BS (k) = ^ J d 3 P P s (p)P B (\\k-p\\), 

which is the Fourier transform of R^(x)Rs(x). From the 
definition of P/?* and Wij(k) we have, 

/>oo 

i f°° 

+ 3%y qP B s{q)dq. 
All that remains is to evaluate, 

Jk J 4tt J Vk(p) 47r||r + p|| 

where Vfc(p) is defined by ||r + p|| > k. We approximate 
this integral using ||r+p|| wpinO<r<p and ||r + p|| ps 
r in < p < r, and find 

/•oo /*oo 
/ qP BS (q)dq « (<5B 2 ) / 9 P S ( g )d g 

/•OO 

+7^ / <?P B (<?)dg. 



Defining, 

IF, 

we have for A 



b/n(*0 = / ^Pb/nO?), 



2I 3 B?m 2 2P rag P- cs w 2 



err ct± 



■ cos 2A 



8M e 2 ff 



P 



where P 2 ,- 



5 l} (SB 2 ) /3. It follows that, 



(A 



2P c ff^ 
2 V M eff m 2 ff 



r 3 



1 / 2P rC gCj 



4 V^effm 2 ff 



cos 



(2A) 



24M 2 ff 



where P 2 ff = P r 2 cg /2 + (SB 2 ) /3. The term that scales as 
L is associated with electron density fluctuations and is 
usually dominant when |A| 3> 1, as is the case for CMB 
photons in galaxy clusters. 

We have assumed that near k = fc cr it and for all large 
k, both P B (k) and Pn(&) are decreasing faster than k~ 3 . 
We have 



k- 1 



•^ 2 - 4xl ° PC V100GH Z 



where ^ = uj/2ir is the frequency of the electromagnetic 
radiation. For galaxy clusters n c ~ 10~ 3 -10~ 2 cm -3 
and for CMB photons z/ ~ 30 - 300 GHz. The typ- 
ical critical length scale is then much smaller than a 
parsec, fc~ it ps 10~ 3 -0.1pc. The form and magnitude 
of P B (k) and Pn(&) in galaxy clusters for such small 
scales are not known empirically. We assume that near 
k = fccrit, P B (k) oc Pjsr(fc) oc fc~ n / 3 corresponding to 
three-dimensional Kolmogorov turbulence. We define 
normalization constants Ck and Cn, 



k 2 P B {k) 



2C K ( — 

ft 



-5/3 



fc 2 P N (fc) = 2(2tt) 1/3 C7^- 5/3 , 
for consistency with Refs. [T(J El], where k — lkpc -1 . 



It follows that, 



W B (k c 



W^(k CTit ) 



6C K (k 



-5/3 



6(2tt) 



1/3 



r 2 /V~ 5/3 



Order of magnitude estimates for the normalisation con- 
stants were calculated in (lQ| |: 



2k 5 /3 C K pa (0.27-0.45) (SB 2 ) L 



-2/3 



2(2^) 1 / 3 C 2 n- 2 « (0.27-0.45) ^pi B 2/3 , 
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where (Sn^\ /n^ = In — 1, and Ln is approximated by 

Thus we find that the dominant contribution to the 
photon- ALP conversion rate is, 



TV, « -(0.27-0.45) 



/ _ x -5/3 



^ 5 / 3 L B 2/3 M- ff 2 L 



., (In - 1) . / , r>2 \ ( 2/ n - 1) 



<<>B 2 ) 



(A3) 



Appendix B: CSZ Contribution to the WMAP 
Signal Extracted from the V + W Bands 

In this appendix, we explain our simplistic approach to 
the processing technique of Komatsu et al. in combining 
the raw data from the WMAP V and W bands in order to 
remove contamination from the CMB anisotropies [16j . 

Considering only the contribution from the thermal SZ 
effect and CMB anisotropies, the observed temperature 
decrement around the cluster is expected to be, 

AT{6, v) = AT CMB (6) + AT SZ {6, v), 

where the fluctuations due to CMB anisotropies are fre- 
quency independent. The frequency dependence of the 



thermal SZ effect was given in Eq. ©. Subtracting the 
signal in the W band from that in the V band we would 
expect (allowing only for the thermal SZ effect) that, 



The results presented in [l6| are expected to be of 
ATl+ w = AT sz , RJ (0). 

However if we include the contribution from the CSZ 
effect, then we expect the total temperature fluctuation 
in any frequency band to be given by, 

AT{6, v) = AT CMB (0) + AT SZ (9, v) + AT CSZ (9, v), 

where ATqsz was given in Eq. ([9]). This leads to a 
predicted signal of, 

ATl+ w (9) = AT sz , RJ (0) 

+AT™ L {6) ( ~ 2 (gcsz ^ gggz) 



5csz(^o)(3sz ~9sz) 



(Bl) 



where g C sz{v) = u 5 / 3 (e~ x - l)/x and v = 204 GHz. 
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